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Abstract 

The aim of this paper is to perform a deeper geometric analysis of 
problems appearing in dynamics of affinely rigid bodies. First of all we 
present a geometric interpretation of the polar and two-polar decomposi- 
tion of affine motion. Later on some additional constraints imposed on the 
affine motion are reviewed, both holonomic and non-holonomic. In par- 
ticular, we concentrate on certain natural non-holonomic models of the 
rotation-less motion. We discuss both the usual d'Alembert model and 
the vakonomic dynamics. The resulting equations are quite different. It is 
not yet clear which model is practically better. In any case they both are 
different from the holonomic constraints defining the rotation-less motion 
as a time-dependent family of symmetric matrices of placements. The 
latter model seems to be non-geometric and non-physical. Nevertheless, 
there are certain relationships between our non-holonomic models and the 
polar decomposition. 

Keywords: affine motion, polar and two-polar decompositions, Green 
and Cauchy deformation tensors, non-holonomic constraints, dynamical 
symmetries, d'Alembert and Lusternik variational principles, vakonomic 
constraints. 



1 Affine constraints, geometry of the polar and 
two-polar decompositions 

Let us begin with a short review of our earlier results concerning the mechanics 
of affinely- rigid body [23, 24, 25j. To be honest, some of them are also partially 
contained in Eringen's theory of micromorphic media, i.e., continua of infinites- 
imal affine bodies [5]. Later on, we developed the theory in various aspects 
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and some of our results were confirmed and developed by many people |15[ 1161 
[171 |2T1 123 SOI SH Let us also mention the papers like [3 H 03 EH OS . 
Nevertheless, in spite of numerous applications the topic does not belong to 
commonly known standards, and because of this a brief repetition seems to be 
necessary. 

Let us consider a system of material points moving in n-dimensional physical 
space M; we assume M to be an affine space with the linear space of translations 
V , endowed also with the symmetric and positively-definite metric tensor g G 
V* ® V*. If necessary, the translation vector from x G M to y G M will be 
denoted by xy. The material space, i.e., the set of material points will be also 
an affine space A of the same dimension n, with the linear space of translations 
U . The material metric tensor will be denoted by r\ G U* <£> U*, and translations 
vectors by at for a, b G U. As usual, we say that a mapping tf> : N — > M is 
affine if it preserves all affine relationships, i.e., there exists a linear mapping 
L [<j>\ : U — > V, denoted also as Defy G L (U, V) such that 

<t> (a) cj)(b)=L [0] at (1.1) 

for any pair of material points, a, b G N. If y % . a K are affine coordinates 

respectively in M and N , this means obviously that <fi is analytically given 
by first-order polynomials: 

y l = x l + ^ K a K . (1.2) 

Obviously, this definition is valid for any, not necessarily equal dimensions 
of N, M . The set of all affine mappings of N onto M will be denoted by 
Af f (N, M), and the set of all one-to-one affine mappings of N onto M is 
denoted by AffI(N,M) (affine isomorphisms). Obviously, AffI(N,M) is 
non-empty only if dim N = dimM, and for any <f> G Af fl (N, M), ip = L [(/>] G 
LI(U,V), i.e., it is a linear isomorphism of U onto V. The groups of affine 
and linear isomorphisms of M and V will be denoted by GAff (M), GL(V). 
They are open subsets of Af f (M), L (V), i.e., of the sets of all affine and linear 
mappings of M and V into themselves. 

Every choice of affine coordinates a K , y l in N, M pre-assumes two things: 
a choice the origins O G N, o G M of coordinates in N, M and a choice of 
bases (. . . , Ea, ■ ■ ■), (• • ■ , &%, ■ ■ •) in U, V, or equivalently, a choice of dual bases 
(. ..,E A ,...), (. ..,e\...) in U*, V*. Then we have 

a K (P) = (E K ,Dp) , y>(p) = (e\4) (1-3) 

for any points P G N, p G M. When the constant co-moving mass distribution 
in N is fixed and described by positive measure fi on N, then it is natural to 
choose D G as the centre of mass, 

J Dpd f i = 0. (1.4) 
2 



The point D is uniquely defined when m — fi(N) is finite, what is physically 
always assumed. With such a choice of O, the quantities x % in (|1.2[) are the 
current coordinates of the centre of mass in M, — (j>{D). Let us stress that 
for any, not necessarily affine, configuration is defined by the condition 

o^pdfi^p) = 0, (1.5) 

where denotes the ^-transport of the measure fj, from N to M. The equality 
= (j) (£)) holds only for affine configurations. 

When the choice of D is fixed as above, then the configuration space of 
affinely- rigid body, i.e., the manifold of affine isomorphisms of N onto M, 
AffI (N, M) becomes canonically identified with the Cartesian product M x 

Li(u,vy. 

<t>={<t>{Q),L[<t>]) = (...,x\...;...,<p % K ,...). (1.6) 

This is the splitting of degrees of freedom into translational and internal ones. 

All those concepts are purely affine and the metric tensors g, r\ occur only 
on the dynamical level. Let us mention, there are also purely affine, metric-free 
dynamical models |34[ 135] , but it is quite a different story. 

The affine groups GAff (M), GAff (N) act on the configuration space of 
affine body AffI (TV, M) through the left and right superpositions. Namely, any 
(A, B) e GAff (M) x GAff (AT) transforms the configuration (f> e AffI (N, M) 
as follows: 

<f)^Ao<f)oB. (1.7) 

Obviously, the action of GAff (M) does commute with that of GAff (N). They 
are respectively the spatial and material transformation groups. 

Describing affine configurations as in (|1.6[) we can represent the actions 71) 
of A e GAff(M), B E GAff(N) as follows: 

A S GAff(M) : (x,(p)->{A(x),L[A]<p), (1.8) 



BeGAff(N): {x,ip)^{t[(p-DB(D)\,tpL[B\), (1.9) 

where for any v £ V, u G U , the symbols t[v], t[u] denote translation operations 
in M and N, i.e., such affine transformations of M and N that 

a;t[w](xj = v, at[u](a) = u (1-10) 

for any x € M, a G N. If, after the material origin D G N is fixed, B is identified 
with (B,b) G GL(U) x U (semi-direct product), then (|1.9|) becomes 

(B,b): (x,<p)->(t[<pb](x),tpB). (1.11) 
Analytically, (|1.8p and (|1.9p / (II. lip are respectively given by 

(. . . , x\ . ...)->(..., A\ n x m + a\ . A? m <p m K , . . .) , (1.12) 
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( -•'....:...,;'„....) • ( V • tp i M x M , /'' k • ...). (1.13) 

The structural difference between spatial (Eulerian) and material (Lagrangian) 
transformations is easily seen here. Let us observe that GL(V), GL{U) act also 
on the manifold Q int = LI(U, V) of internal/relative degrees of freedom through 
the obvious formulas: 

tp -> ApB , (A, B) G GL{V) x GL{U). (1.14) 

Obviously, this action is non-effective and corresponding kernel is given by the 
subgroup: 

{(Xldv^^Idu) : AgR\{0}} C GL(V) x GL(U). (1.15) 

The subgroups of GL(V), GL(U) and those of GAff(M), GAff(N) act in a 
natural way on the internal and total configuration spaces Qi n t — LI(U,V), 
Q = Aff(N,M) ~ M x Qi n t- Let us mention a few most important of them: 
orthogonal groups 0(V, g), 0(U, rj) their rotation subgroups SO(V, g), SO(U, rj), 
special linear groups SL(V), SL(U) or one-dimensional dilatation subgroups 
Dil(V) = {XId v :AeR \ {0}}, Dil(U) = {XId w : A G R \ {0}}. In the total 
configuration space Q, when translational degrees of freedom are taken into 
account, those groups are semi-directly extended by translations T(M) ~ V, 
T(N) ~ U to the corresponding affine subgroups: Euclidean E (M, g), E (N, 77), 
isochoric SAff(M), SAff(N) and dilatations/translations Dil(M), Dil(N). 
The meaning of symbols is obvious. Let us only remind a few definitions. A G 
0(V,g), B G 0{U,rf), tp G 0(U,ri;V, g) when they preserve the metric tensors, 
thus, 

g%j = 9klA k iA l j, t/ab = i] C dB C aB D b, Vab = gij^A^B- (1-16) 

A e SO(V,g), B £ SO(U,rj), when not only |det A| = |detS| = 1, but just 
detA = deti? = 1. When orientations p, uj in U, V are fixed and dettp = 1 
in some orthonormal positively oriented bases in U, V, then we say that ip G 
SO (U, 77, p; V, g, uj) when det tp = 1 and that tp e O (U, 77; V, (7) when |det y| = 1 
Similarly, we say that A G SL(V), B G 5L(?7) when detA = det 5 = 1 
but without orthogonality condition (|1.16[) . And similarly tp G SL {U, p; V, ui) 
when det tp = 1 in some positively oriented bases and that tp. G UL(U, V) when 
I det tp\ = 1. If I det A\ = |det B\ = 1 we say that, A, are unimodular and write 
that A G UL(V), B G UL(U). If det A = det i? = 1 we say that they are special 
linear. 

Lie algebras of GL + (V), GL + (U), the proper (positive-determinants) sub- 
groups of GL(V), GL(U) are isomorphic with the commutator Lie algebras 
of all linear mappings GL(V)' ~ L(V), GL{U)' ~ L(U). And Lie algebras 
SO(V, g)', SO{Urj)' consist respectively of g- and 77-skew-symmetric elements of 
L(V), L(U): 

4 = -a/ = -g jk g il a\ b A B = -b B A = - V BcV AD b c D . (1.17) 
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Lie algebras SL(V)' , SL(U)' consist of trace-less linear mappings: 



Tra = a\ = Q, Trb = b K K = 0. (1.18) 

As usual, when dealing with group-theoretic degrees of freedom, it is convenient 
to use non-holonomic Lie-algebraic velocities. In the case of affine systems we 
shall use the term "affine velocity"; Eringen referred to them as "gyrations". The 
spatial and material affine velocities of internal motion are given by 

n = i- 1 6 l(V), n = cp- 1 ^ = ^nip e l(u). (i.i9) 

at at 

Besides the usual velocity of translational motion, v % = dx 1 /dt, one uses also its 
co-moving representation: 

d=ip- 1 v, d A = (ip- 1 ) A l v l . (1.20) 

When gyroscopic constraints of metrically-rigid motion are imposed, 

<peSO{U,t],p-,V,g,u), (1.21) 

then fl, f2 are respectively g- and 77-skew-symmetric, i.e., they satisfy (|1.17p when 
substituted instead a, b. This is the alternative, "anholonomic" representation 
of those holonomic constraints. 

When the body is incompressible, dettp = 1, then f2, Q are trace-less, i.e., 
they satisfy (fO| . 

Gyroscopic and isochoric constraints are, obviously, holonomic. Geomet- 
rically this has to do with the fact that their affine velocities are elements 
of the commutator Lie subalgebras of L(V), L(U). There are, however, an- 
other interesting cases of non-holonomic constraints of spatially and materially 
rotation-less motion. In the first case f2 is g-symmetric, in the second one Vt, is 
^-symmetric, i.e., respectively, 

n* 3 - = n/ = g jk g a n k h n A B = n B A = riBcv AD ^ c D. (1.22) 

Let us observe that unlike in the holonomic gyroscopic constraints, the two 
conditions (|1.22p are non-equivalent and describe different non-holonomic con- 
straints. Namely, the first, i.e., spatial, condition in (|1.22ll is materially repre- 
sented by 

& A b = G BC G AD h c Dl (1.23) 
where Gkl are components the Green deformation tensor, 

G K L= Oij^K^L, (1-24) 
and G KL represent its contravariant inverse, 

n KM r _ cK r KL 1 KM LN /~i /-, 9 r\ 

G LrML — "Li G f= rj r\ Unm- U-^oJ 
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Similarly, the second, i.e., material, equation in (|1.22p is in the spatial language 
equivalent to 

f2j = C Cj m n m k, (1.26) 
where Cy are components of the Cauchy deformation tensor, 

C id =VABW- l ) A i{v- 1 ) B h (1-27) 

and C % i are coordinates of its contravariant inverse: 

C ik C kj = 6% C ij y£ g lk g ]l C kl . (1.28) 

The inequalities in (|1.25l) . (|1.28p show that the thoughtless use of the kernel- 
index convention may be misleading. Having in U, V two metric-like tensors G, 
rj and C, g one can construct two mixed tensors: 

G A B := V AC G CB , &,j := g ik C kj (1.29) 
and the family of scalars, e.g., 

h = Tr (G fc ) = Tr (cT fc ) . (1.30) 

Those scalars are invariant under the action (I1.14[) of the subgroup 0{V, g) x 
0(U,rf) C GL(V) x GL(U). They are basic orthogonal deformation invariants 
of ip. According to the Cayley-Hamilton theorem, there are only n independent 
invariants, e.g., those corresponding to k = 1, . . . , n\ any other invariant is their 
function. 

Deformation invariants tell us how strongly the body is stretched/contracted, 
but they do not contain any information about the spatial or material orientation 
of the stretching. This information is encoded in directions of the main axes 
of deformation tensors G, C. More precisely, let L a [<£>], R a [<f], a = 1,... ,n, be 
orthonormal basic eigenvectors of C[ip], G[ip]: 

CL a ^X a 1 L a} GR a = X a R a , (1-31) 
g {L a , L b ) = gijV'aUb = $ab = VcDR C a R D b = V {R a , Rb) ■ (1-32) 

Their dual covectors L a [ip\ e V* , R a [ip] £ U* satisfy 

CM = X>a X M£>]®£ a M, (i.33) 

a 

GM = J2 X °M Ra M® Ra M- (1-34) 

a 

It is convenient to introduce the symbols Q a , q a , 

Q a = cxp (q a ) = yfc. (1.35) 

The ordered bases £=(..., L a , . . . ), R = {. . . , R a , . . . ) are identified with 
isomorphisms L : R" — > V, R ; K™ — > U and their dual co-bases L" 1 = 
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(. . . , L a , . . . ), R^ 1 = (. . . , R a , . . . ) may be interpreted as the inverse isomor- 
phisms L _1 : V — > R™, R^ 1 : U — > M". The diagonal matrix with diagonal 
entries Q a , Diag (. . . , Q a , . . . ), is identified with an isomorphism D : E™ — > W l . 
And finally, the isomorphism ip : U — > V may be represented as: 

tp = LDR- 1 . (1.36) 

In this way affine configurations ip are identified with the triplets consisting 
of two gyroscopic configurations of metrically-rigid bodies L, R, and of the 
system of n material points on R (deformation invariants). 

For any pair of linear bases, e.g., (. . . , i? a [</?], . . . ) in U and {. . . , L a [ip], . . .) 
in V, there exists exactly one linear mappings U[(p] of U onto V such that 

L a [<p]=U[<p]Ra[(p], a = l,...,n. (1.37) 

When the bases are respectively rj- and g-orthonormal, then U[(p] £ 0(U, rj; V, g), 
and ip may be expressed as 

<p = U[<p]A[<p] = B[<p]U[<p], (1.38) 

where the linear mappings A[ip] € GL(U), B[tp] G GL(V) are symmetric in the 
rj- and <?-sense and positively-definite. Obviously, by the positive definiteness 
we mean 

77 (A [ip] z, z) > 0, g(B[<p]uj,uj)>0 (1.39) 
and by the symmetry 

v(A[(p]u,v) = tj(u, A [(p]v) , g (B [<p] x, y) = g (x, B [p] y) (1.40) 

for non- vanishing z € U, u> € V and for any it,u S C/, cc, y S V. 

Unlike the two-polar splitting (|1.36p the both versions of the polar splitting 
(|1.38p are unique. And, obviously, A[ip], B[(p] are related to each other by the 
[(^-similarity : 

A[<p]=U[<p]- 1 B[<p]U[<p]. (1.41) 

Analytically, in the matrix language L, R, U are orthogonal, D is diagonal 
positive, and A, B are symmetric positively definite matrices. 

Transformations (|1.7[) , (|1.14|) act on affine velocities according to the obvious 
rules: 

n^AnA' 1 , h^B^hB. (1.42) 

According to the same rules orthogonal transformations act on the g- and rj- 
skew-symmetric angular velocities. 

In analogy to non-holonomic affine velocities, one introduces their dual affine 
spin quantities, 

E = tpP, E = Pcp = ^ _1 E<^, (1.43) 

where P denotes the system of canonical momenta P\ conjugate to ip % A, E^, 
Tj A b are momentum mappings of the transformation group (|1.14|) . In other 
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words, they are Hamiltonian generators of this group. Their Poisson brackets 
correspond to the structure constants of GL(V), GL(U): 

{E^E^} = Sh&j - S^TS, (1.44) 
[t A B ,YF D ] = 8 c b % a d -8 a d YF b , (1.45) 

{e^S^b} = 0. (1.46) 
And, obviously, the following holds: 

A} = SWa, {x A B ,<P k c} = 8 A c<p k B . (1.47) 
Transformations ([TT]) . (fl~T4"f act on S^-, S A B just like on Q, Q ([T32]) : 

E — > AEA -1 , %^B- 1 %B. (1.48) 

The same may be done with the gyroscopic degrees of freedom of the two-polar 
and polar decompositions. The R"-comoving angular velocity x% of the L-top 
and the F-spatial representation \ l j are given by 

M^H^p <L49) 

AT i 

X i j = ^L a j , X = x\L a ®L b . (1.50) 

And similarly, the R n -comoving and [/-spatial components of the angular ve- 
locity of the i?-top, {) l are given by 

A = d = d\R a ®R b . (1.52) 

ere 



Obviously, the "comoving" angular velocities x% are (5-antisymmetric, while 
the "spatial" ones, x i ji $ £> are respectively g- and 77-antisymmetric. The dual 
spins conjugate to x%, Xji are also skew-symmetric matrices, denoted 
respectively by j5°b, t%, p l j, t l- When canonical momenta conjugate to q a are 
denoted by p a , then the duality relations have the form: 

((p,T,p), \X,$,cn} = ((p,r,p),(x, 

= Pa<f + T^Tr (px) + ^Tr (ti?) . (1.53) 
It is clear that p, r are Hamiltonian generators of the transformation groups, 
ip^AtpB-- 1 , AeSO{V,g), B e SO(U, 77), (1.54) 
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so, they are equal respectively to the spin and minus vorticity. Let us remind 
that the spin and vorticity are doubled g- and 77-skew-symmetric parts of £*,-, 
Tj A b- When we use the polar splitting (|1.38[) . then the gyroscopic /7-motion is 
characterized by angular velocity in the co-moving and spatial representations, 
respectively: 

Q = U-^, u^^U- 1 = UQU-\ (1.55) 
dt dt v ; 

2 Kinetic energy, equations of motion, additional 
constraints 

It may be easily shown that the kinetic energy of the classical affine body is 
given by 

T = T tr + T mt = ^ 9ij d 4 d 4- + Xm^^T-J*** (2-1) 

where m G K and J G U <8> U are constant inertial parameters of affine degrees 
of freedom, 

m = J dfi(a), J AB = J a A a B dfi(a). (2.2) 

N N 

Therefore, m is the total mass of the body and J AB is the quadrupole momen- 
tum of the mass distribution, algebraically equivalent to the co-moving inertial 
tensor. Let us repeat that, following (|1.4j) . the dipole momentum vanishes: 

a A dfi(a) = 0. (2.3) 

And the higher multipoles, although no 11- vanishing, do not contribute to the 
affine motion. 

The kinetic energy (I2.1[) is invariant under the action (|1.14l) of 0(V,g) x 
0(U, J -1 ) C GL(V) x GL(U), and of course under the action of translations. 
It is important that because of the essential dependence of its metric tensor on 
g, J, it fails to be invariant under the total action of affine groups GAff(M), 
GL(V), GAff(N), GL(U) .After substitution of the polar decomposition (fi~38)) 
to (]2.1|) . it becomes 

1 dA K A dA L B AB K c dA L B AB 

lint = T^VKL — T, -J— J cA A J 

2 at dt dt 

+ l -mLU K c O L D A c A A D B J AB . (2.4) 

The first term represents the kinetic energy of deformative vibrations, the second 
one is the Coriolis coupling between deformative and rotational motion, and 
the third term describes the centrifugal coupling of rotations and deformations. 
The lower-case indices of the third term are contracted with the A-deformed 
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inertial tensor, A c aA d b J ab ■ A similar formula holds for the second of the 
polar decompositions (|1.38p . When we use the purely analytical language and 
orthonormal coordinates, t)kl = 5k l, then the formula (|2.4[) may be written in 
the following brief matrix form: 



T - = \ T - ( J ) J + Tr { A/ i Q ) - ¥ < AJAQ "> ■ <"» 

Substituting the two-polar decomposition (|1.36l) and its by-products (|1.49[) . 
(|1.50p . (|1 .51)1 . (|1.52p to (|2.ip , we obtain a rather complicated formula: 



^Tr (D^DR^JR) - ^Tr (dD 2 dR~ 1 JR 
1„ A^„~„~_ 1td \ !™ / n-n.4D-l 



+ -Tr [dDxDR-'JRj + -Tr [DxD'&R~ 1 JRj . (2.6) 

It is seen that the complication is due to the term R~ X JR. And indeed, the 
to polar splitting is computationally optimal in the special case of inertially 
isotropic body, when 

jAB = Irj AB (2 ?) 

and the kinetic energy is invariant under the material orthogonal group 0(i7, rf). 
Then 

(R- 1 JR) ab = IS ab (2.8) 
and the formula (|2.6p simplifies to 



T„ t = -Tr ^-j J + ITr \DxD&) - -Tr (D 2 x 2 ) - -Tr [D 2 § 2 ) . (2.9) 

It is clear that combining appropriately x> $ we can avoid interference terms. 
This is explicitly seen when instead of the usual kinetic formulas like (|2.ip , (|2.4p , 
(12. 5p , (|2.9p one uses their canonical forms based on the Legendre transformation, 
like, e.g., 



dx J A dip- 



Pi = m 9ij =^, p A L ^ gi] ^J AB (2.10) 

in the case of the usual, velocity-independent potentials. Then instead of (|2.ip 
we obtain 

T =T tr + T int = ^hiPi + ~ (J- l ) AB v\v B i9 l \ (2.11) 
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where, obviously, g %: > is the contravariant inverse of gtj, and (J x ) Ag is the 
covariant inverse of J AB , 

//•'.<//.,. <>V (J- 1 ) AC J CS =^ S . (2.12) 

Obviously, if there is a velocity-dependence in the potential, then the formula for 
T is more complicated. For example, the presence of magnetic fields results in 
the configuration-dependent translational gauging of canonical momenta. But 
in a moment we are not interested in such details. Making use of the duality 
(| 1 . 53[) we can write the Hamiltonian form of (|2.9p as follows: 

T 1 2 1 (M%) 2 1 ^ ( Na bf 

mt " ^V a + ^hW^f + ^(^TW* (2 - 13) 

where the quantities M%, N% are given by 

M% := -p?b - r\ N% := jft - r% (2.14) 

Q a are given by (|1.35p . and P a are their conjugate momenta, 

P a = p a cxp (-q a ) . (2.15) 

It is seen that in (|2.13[) one deals with a kind of "diagonalization" of the expres- 
sion for Tint ■ 

Equations of affine motion may be derived on the basis of the variational 
principle for the following Lagrangian: 

L = T — V(x, ip), (2.16) 

where T is given by (12. ip . or in the Hamiltonian terms: 

^ = {F,H}, (2.17) 

where H is the Hamiltonian corresponding to L, and F runs over a set of 2n(n + 
1) independent phase-space functions. But independently of this variational 
framework, they may be derived in general, on the basis of d'Alembert principle. 
According to this principle, equations of affine motion are obtained from the 
general equation of motion of the underlying system of material points by taking 
the monopole and dipole moments of the balance laws for the linear momentum. 
The points is that the original equations should be modified by introducing the 
reactions responsible for maintaining of the constraints. The reactions them- 
selves do not vanish, but their monopole and dipole moments do so. Because of 
this, the effective, free of unspecified reactions equations of affine motion have 
the form of the balance laws for the total linear momentum and the total affine 
momentum (hypermomentum, affine spin): 

dtf =pi dK^ = d^d^ jAB + NlJ 
dt dt dt dt 
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with the following meaning of symbols: 

k = J v{P)dfi(P), (2.19) 

K = j O^P)^) ® v(P)dfi(P), (2.20) 

N = J o4P)<P(p] ®F(P)dn(P). (2.21) 

Therefore, k is the total linear momentum, K is the total dipole moment of 
the momentum distribution (taken with respect to the instantaneous position 
of the centre of mass) and N is the total dipole momentum (affine torque), also 
related to the centre of mass instantaneous position. Reaction forces responsible 
for affine constraints are automatically cancelled in (|2.19l) - (|2.2ip . The formulas 
(|2.19[) — (|2.21[) are general, but their affine versions are just 

dt r dt 

in the sense of symbols (|1.2[) . 

Let us observe that the usual Legendre transformation (for velocity-indepen- 
dent, usual potentials) identifies those kinematical quantities with Hamiltonian 
ones pi, E',- up to the index position, inessential in Cartesian coordinates: 

Pi </,>'. E*,- =K im g mj . (2.23) 

It also convenient to use the co-moving representation of linear momentum, 
affine spin and affine moment of forces: 

(2.24) 

l ) B .K ij , (2.25) 

- 1 ) 3 N* 3 . (2.26) 

As it was said, the general, non-variational equations of motion are given by 
(|2.18p . (I2.22p . Let us quote a few equivalent forms also based on (|2 . 18[) with 
substituted (|2.22l) . like, e.g., 

dk l dK ij 

-7T = F\ n\ n K^ + N" . (2.27) 

dt dt 

Let us also observe that one can write: 

^_ =N t 3+2 ^HL. (2.28) 
dt ogtj 

In particular, for the doubled skew-symmetric part of K 13 , i.e., for the angular 
momentum, we have 

dS' 1 dK' i dl\ > ' 

^— = ^—-^—=N l 0-N^=Af'\ (2.29) 
dt dt dt 



P A = 






k AB = 






N AB = 







12 



where the right-hand side denotes the usual torque. If A™ is symmetric, in 
particular vanishing, this becomes the usual conservation of angular momen- 
tum. The purely Lagrangian, i.e., [/-based form of equations of motion, may be 
formulated as follows: 



dk A 
~di 
dK AB 



- -k B (J- 1 ) BC K CA +F A , (2.30) 

= -K AC (J- 1 ) CD K DB + N AB . (2.31) 
Expressing these equations in terms of kinematical quantities one obtains 

m— = -mn A B v B +F A , (2.32) 
dt 

^^J CA = -n B D n D G J CA + N AB . (2.33) 
The explicit form of equations of motion reads: 



dt 2 V ' * dt 

ipA ^ 2 ~ J = N [ x '-a** '-*-'*)- (2 - 35) 

The assumed non-singularity of matrices W> l A \ , [J AB ] m principle enables 
one to solve the second equation with respect to second derivatives d 2 f l A/dt 2 , 
expressing them through dynamical variables. Nevertheless, the form (|2.35|l is 
more convenient, because it is geometrically suited to the nature of our prob- 
lem, in particular to additional constraints which may be imposed on the afhne 
motion. This follows from the fact that in affinely-rigid behaviours the formula 
for the power of forces is given by 

V = V tr + Vint = Fig l3 v* + \ !! !i,iM[, = FA + Ni^ ■ (2.36) 

And if for F J , N 3 one substitutes reactions maintaining constraints, this ex- 
pression vanishes; only the external given forces contribute here. Similarly, if we 
subject the general affine motion to some natural group-theoretical constraints, 
in a consequence of which, e.g., O does belong to some Lie subalgebra of L(V) 
or to some other linear subspace of clear algebraic meaning, then the effective 
reaction-free system of equations of motion consists of some natural subspace 
of (|2.35p and, obviously, of the explicit description of constraints. This would 
not be the case if we used the form of ()2.35[) solved with respect to the second 
derivatives of ip>B- Let us quote a few convincing examples of both holonomic 
and non-holonomic constraints. 

2.1 Metrically rigid motion 

It consists in that both the mappings 0, ip are metrical isometries, therefore, 
all distances and angles are preserved during the motion, so that the following 
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holds: 

Vab = gijf 1 a^b- (2.37) 

This explicit holonomic representation implies that fi i j , £L A b are respectively g- 
and fy-skew-symmetric during any admissible motion, 

Wj = -Slf = -9ja9 ib n%, tt A B = -&b A = -VBCV AD & C D- (2.38) 

Therefore, (|2.35|) is not valid any longer, instead the following equations with 
unspecified reaction torques hold: 

= -^JT-^J + N 3 + Nr 13 , 2.39 
dt dt dt 

or equivalently 

^a^-J AB = + N R ij . (2.40) 

But the d'Alembert principle tells us that the power Vr of gyroscopic reactions 
vanishes on every Q l j compatible with constraints, i.e., g-skew-symmetric: 

Vr = Nn^tfi = N R ij n J( = (2.41) 

if (|2~55f holds. Therefore, 

N R ij = Nr* , Nr^ = N R / = g ja g ib N R \. (2.42) 

But this simply means that the reaction-free gyroscopic equations of motion 
consist of the skew-symmetric part of (|2.39|) or (|2.40[) with substituted (|2.37p . 
e.g., 

tfJ^BjAB _ ^fgBjAB = N ij _ N ji = M i3_ (2.43) 

dt dt z 
In other words, any parametrization of the isometry manifold ("Euler angles", 
"rotation vectors", etc.) may be safely substituted to (|2.43p . The torque N 1 -* = 
N l ° — N^ 1 becomes the function of those parameters and we obtain the system of 
n (n — 1) /2 independent of equations (j2.43[) imposed onn(n- 1) /2 parameters. 
It is clear that this system is the balance law for spin: 

d f*J d 

— = — (K ij - K ]l ) = Af lj . (2.44) 
dt dt y ' K ' 

It becomes the spin conservation when the torque J\f lJ does vanish, i.e., when 
jV ij is symmetric. One can as well rewrite (|2.44p in Lagrangian terms: 

is. 



(k AB -K BA ) = K BC (J- 1 ) ( 



dt dt\ ' y > CD 

^CD 

where 



- K AC {J- 1 ) CD K DB +M AB , (2.45) 



N AB = N AB ~N BA . (2.46) 
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This may be easily expressed in terms of co-moving affine velocity: 

d^CjCA _ dn^C jCB = ^A d ^D cJ CB _ ^B^D^CA + ftAB {2AJ) 

Those are Euler equations. When the co-moving inertial tensor is spherical, i.e., 
when 

J = -7?, (2.48) 

then the non-dynamical terms on the right-hand side of (|2.47[) do vanish, and 
we obtain simply 

dh AB dh A 



, =/Zr_£ ?? C7S = ftAB (2.49) 

2.2 Shape-preserving motion 

Now the shape of the body is preserved, but not necessarily its size, so that 

gij<p l A<P 3 B = ^VAB, (2.50) 

where A denotes the time-dependent coefficient. Then, in analogy to the previ- 
ous example, the d'Alembert principle tells us that the react ions- free equations 
of the constrained motion consist of the skew-symmetric part and the trace of 
or (12351) : 

d 9 d 

= — (K ij ~ K ji ) = N 13 = N ij - N ji , (2.51) 

dt dt v ' K ' 

^ - ±(*#>)-»4£*£j» + »S*-*r + *. CM, 

when written in a few independent forms. Using directly the representation in 
terms of coordinates, we obtain that 

tfJ^LjAB _ v i A f^_jAB = M ij = N ii _ ^ji (2 53) 

dt A at* 

9i ^ A ^w~ jAB = 9ijNij - (2 - 54) 

2.3 Incompressible afRne motion 

Incompressibility (isochoric motion) means that 

- det [ip\] = 0. (2.55) 

This condition is well defined, although det \}P % a\ is not a scalar but scalar 
density with respect to both spatial and material coordinate transformations. 
The identity 

| det [^a] = det [pi B ] {v^ti^f = ( 2 - 56 ) 
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is equivalent to 

Trn = Trn={ip- 1 ) A i ^=Q. (2.57) 

D'Alembert principle implies that reactions Nr which keep these constraints, 
being dual to the subspace of all trace-less matrices, 

V R =N R i jrti = Q, fi** = 0, (2.58) 

are proportional to the identity mapping, 

W*',- = A<$*-, A/," A//'-'. (2.59) 

Therefore, the effective reactions-free system of equations of motion is given by 
the trace- less part of the original balance law for , i.e., explicitly 

i T AB 1 « ^Vb TAB„ij AT ij 1„ AT abij (1 ad\ 

<P A dt2 J ~ n b ~dt 2 ~ = ~ n 5 (2.60) 

This is a system of (n 2 — 1J independent equations of motion imposed on 
(n 2 — l) independent parameters of ip l A (cf. (|2.55I) V 

2.4 Spatially rotation-less motion 

We have seen above that the concept of the purely rotational degrees of free- 
dom is well defined and correctly formulated: simply the mappings tf>, (p are 
isometries. It is not so with the opposite concept of rotation- free, i.e., purely 
deformative configurations. The first, naive idea would be to base this concept 
on the polar decomposition, either in left or right version. So, rotation-free con- 
figurations would be ones given by the purely deformative factor in the polar 
decomposition. Therefore, depending on whether one deals with the [/-left or 
[/-right version of (|1.38| . we would say that tp is purely deformative when it 
coincides with its r\- or ^-symmetric part A[ip] or B[(p]. However, this would be 
incorrect. The more incorrect would be attempts of introducing the pure defor- 
mation on the basis of the two-polar decomposition (|1.36[) . There are a few deep 
geometric reasons for that. First of all, the symmetric mappings A[ip], B[ip] do 
not describe any configurations at all. It is only mappings from U to V that may 
be used as a model of the configuration, neither the linear automorphisms of U 
nor those of V. Without fixing some standard element of the manifold of isome- 
tries O ([/, 77; V, g) we cannot identify automorphisms of the material or physical 
spaces with any mappings of U on to V. So, even from this relatively naive 
point of view, the symmetric mappings do not describe configurations. But 
there are also other arguments. Namely, even if we "forget" about the above 
fact and simply proceed with the R™-model of space and body manifolds, using 
the elements of GL(n,M) as LI(U, V), it is still so that the symmetric matrices 
do not form a Lie group. Therefore, the corresponding absence of rotation is not 
an equivalence relation because of the transitivity failure. If some configurations 
A, B are mutually non- rotated in the polar sense, i.e., they are related by the 
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symmetric matrix Sym(A, B), and if so are B, C in the sense of being obtained 
from each other by the action of some Sym(B, C), then in general A, C are not 
connected by a symmetric matrix. It is well known that the symmetric matrices 
do not form a Lie group or Lie algebra. The product Sym(A, B) Sym(B, C) 
in general is not symmetric; instead it splits into the multiplication of some 
symmetric matrix and some nontrivial isometry. So, certainly, being related by 
a symmetric matrix is not an equivalence relation, and the concept of mutually 
rotation-free configurations is not correct. But there are well-defined rotation- 
less motions. We say that the motion I 3 i 4 ip(t) € LI(U, V) is spatially 
rotation-less when is ^-symmetric: 

Wj - 0/ = - g ]k g U ^ k i = 0. (2.61) 

And similarly, we say that it is materially rotation-less when f2 is ^-symmetric: 

h A B - n B A = n A B - VB cv AD n c D = o. (2.62) 

The symmetry of fl or Q is just the natural complementary concept of their 
antisymmetry in rigid motion. And therefore, this is a proper definition of the 
rotation-less behaviour, just behaviour, not configuration. The point is that 
the symmetric matrices do not form a Lie algebra. On the contrary, they are 
anti-Lie algebras in the sense that their commutators are respectively g- and 
77-skew-symmetric: 

[Sym(L(V),g),Sym(L(V),g)] = Asyra (L (V) , g) ~ SO (V, g)' , (2.63) 
[Sym (L (U) , r,) , Sym (L (U) , rj)} = Asym (L (U) , 77) ~ SO (U, rf}' . (2.64) 

This is an interesting example of non-holonomic constraints, in a sense different 
than the classical constraints of slide-free motion. Nevertheless, some relation- 
ship with the usual non-holonomic problems of non-sliding motion still seems to 
exist in certain hypothetical applications. Let us consider, e.g., an affine motion 
of a small inclusion or droplet suspension in very viscous fluid. It is natural to 
expect that the surface friction may be an obstacle against rotations. And then 
probably the effective constraints of rotation-less motion may appear. 

Let us stress some circumstance. Namely, the holonomic gyroscopic con- 
straints may be written alternatively in two apparently non-holonomic forms: 

n) + n/ = 0, h A B + h B A = o. (2.65) 

They are mutually equivalent. On the other side, the two versions of non- 
holonomic constraints (|2.61[) and (|2.62p are non-equivalent. Namely, the g- 
symmetry of fl is equivalent to the G-symmetry of fi where, as usual, Gab 
denotes the Green deformation tensor, so that (|2.6ip is identical with 

Gac&b - G B c& c a - 0. (2.66) 

In a moment we are unable to answer the question concerning the details of 
this relationship and the possible fields of physical applications. From a perhaps 
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naive point of view, it seems to be so that it is rather the Euler symmetry (|2.6ip 
that seems to be applicable to description of the afhne motion of suspensions in 
viscous fluids. 

In any case, it is an interesting and rather new problem to discuss the struc- 
ture of equations of motion subject to rotation-less non-holonomic constraints. 
Again the d'Alembert principle shows the advantage of the if-balance form 
of equations. Namely, the effective, reactions-free equations are given by the 
symmetric part of the balance laws, 

+ ^ d ~^ jAB = ^ + (2 - 67) 

together with the algebraically substituted constraints (|2.6ip . The right-hand 
side of (|2.67[) depends only on given forces and is free of reactions. The Lagrange 
form of (|2.67p is given by 

d^C jCA+ dn^c jC B = _fiB D fiD cJ CA_fiA D fiB c jCB + ftAB + ftBA^ (2 g g) 

at at 
where is subject to (|2.66[) . 

2.5 Materially rotation-less motion 

It is also non-holonomic and somehow related to the spatially rotation-less sit- 
uation, nevertheless, in our opinion it is a bit less intuitive. Now the material 
gyration is assumed to be ^-symmetric, i.e., (|2.62p is assumed to hold. The 
effective, reaction-free equations of motion may be written as follows: 

dk AC ~ rlK BC ~ ~ ~ 

^—V C B + ^^V C A = N AC V C B + N BC V C A 
at at 

-K AM (J-i) MN K NC V C B - K BM (J-*) MN K NC V C A , (2.69) 

where, obviously, K AB are co-moving components of K l i , and 

k AB = n B c J AC , V A B = G Ac rf B . (2.70) 

These equations are much more complicated than those for the spatially rota- 
tion-less motion. Namely, their non-dynamical terms depend on the Green ten- 
sor, therefore, also on the configuration ip. The Euler form is also complicated: 

<P l A^j^J AB 9bcC c 3 + ^ A f^J AB g bc C a = N* b g bc C c i + N' b g bc C c \ (2.71) 
where 

Cab ,a b ^AB /o *7<"A 

= ip A (f BV (2.72) 

is the inverse Cauchy tensor. 
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3 Dynamical symmetries of affine motion 



The non-holonomic constraints of rotation- less motion, i.e., the above examples 
described in subsections 12.41 and 12. 51 are really exceptional and in a sense sur- 
prising within the realm of constrained affine motion. Let us stress to avoid 
some easy misunderstandings: they are really non-holonomic and have nothing 
to do with apparently suggestive constraints of the type that (p in (| 1 .38[) is sym- 
metric. Moreover, our analysis above shows that such a formulation would be 
inconsistent, just because of fixing some of infinitely possible isometries U[<p]. 
And, let us repeat, the symmetric matrices do not form a Lie group. The sym- 
metry of or leads to certain equations satisfied by U[tp], A[(p], B[ip] in (|1.38l) . 
but these equations are differential, not algebraic ones. The ^-symmetry of fl 
or fy-symmetry of f2 are the only natural counterparts of their antisymmetry in 
rigid motion. And in any case, they are geometrically interesting special cases 
of constraints, worth to be investigated from the very point of view of purely 
analytical mechanics. 

It is interesting to "solve" the constraints equations (|2.61j) . i.e., to "paramet- 
rize" somehow the manifold of non-holonomic constraints. The best candidates 
are suggested by the polar decomposition (|1.38|) . Let us remind that U[<p] is an 
isometry and that A[ip] is ^-symmetric, thus, 

r) A B=9 ij U{<p} i A U[<p\ j B, (3.1) 

rjAcA c B = VbcA c a, "f] Ac ^-A C B = riBc^-A c A , (3.2) 

at at 

and the co-moving angular velocity of the [/-rotator, uj £ O (U, 77)' C L (U) is 
given by (| 1 . 55[) and is, obviously, 77-skew-symmetric: 

r) A cu} C B = -Vbcu C a- (3.3) 

Substituting those conditions to the definition (|1.19|) of the affine velocity Q, we 
obtain after easy calculations the conclusion that 

_ A 1 ( ( A -i\A dA c B dA A c , , n c 



^_^-. )V ^-^(a-t.J. m 

Therefore, the angular velocity of the [/-rotator equals the half of the commu- 
tator of two algebraically independent instantaneous quantities A^ 1 , dA/dt: 



_ 1 

"=2 



A ' dt 



(3.5) 



In any case, this quantity in general does not vanish and this reflects the non- 
holonomic character of our constraints of non-rotational motion. It is something 
different than the constancy of U, i.e., the vanishing of cD. Making use of the 
polar decomposition (|1.38|l and gyroscopic angular velocity ()1.55j) , we can, a bit 
formally, write down the constraints equations (|3. 511 in the following Pfaff form: 



U~ 1 dU - -A~ 1 dA + - (dA) A^ 1 = 0. (3.6) 
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This system of Pfaff equations is evidently non-integrable. 

Before going any further with the analysis of the constrained affine motion, 
let us quote a few remarks concerning the invariance problems. We are interested 
mainly in the internal, i.e., relative, motion and concentrate on the spatial and 
material rotational invariance. The second, i.e., internal, equation (|2.35p is 
0(V, g)-invariant if for any its solution t tp(t) and for any A s 0(V, g) the 
motion t — > Aip(t) is also a solution. This implies that 

Wa^WuW^SV (3-7) 



dt J V dt , 

where we do not indicate explicitly the possible explicit time-dependence of AT. 
But (|3.7p means that the co-moving representation of N is non-sensitive with 
respect to the action of A £ 0(V, g): 

This means that N is algebraically built of the co- moving quantities G, and 
any fixed material tensor K in U: 

N(v,^)=F(G,fi,K), (3.9) 

It is interesting that this form of N implies the rotational invariance of equations 
of motion, but it does not imply the conservation of spin, i.e., internal angular 
momentum. Spin is conserved only if N is a symmetric tensor, just like in 
mechanics of micropolar or micromorphic continua. Let us mention in particular 
that the Green-Ostrogradskij theorem implies that in orthonormal Cartesian 
coordinates the afhne momentum of forces is proportional to the mean value of 
Cauchy stress tensor in the medium: 



j a l K (3.10) 



In the case if hyperelastic bodies, both continuous and discrete, A™ , N AB are 
automatically symmetric. Indeed, the condition for the potential energy 

V(A<p) = V(<p), AeO{V,g), (3.11) 

implies that V is algebraically built of the Green deformation tensor: 

V{y)=w(G[ip],K), (3.12) 

where K again denotes any state-independent tensor in U. And then one can 
show immediately that 

N( V ) AB = 2j^ = N{v) BA , N«=W\ (3.13) 
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therefore, spin is a conserved quantity. 

This was about the invariance under the left-hand side action of O (V, g) 
on internal/relative degrees of freedom. Let us now ask what are conditions of 
the invariance under the right-hand side action of material orthonormal group 
O (U,r]). One can easily show that for any solution t — > ip(t) of (|2.35p and for 
any B £ 0(U, rj) the right-rotated motion t — > f{t)B is a solution too when the 
following holds: 

d V'L dK oL jCD _ f..D d( P 



^ K J$r L ] cB L D J^ u = N» UpB, -£B) . (3.14) 

Unlike in the case of spatial isotropy, this implies two conditions — the 
internal and dynamical ones: 

J = I V , N(vB,%B)=N( V ,%). (3.15) 

The second conditions in (|3.15|) implies that N depends on the mechanical state 
(ip, dtp/dt) through the pair (C, 51) and any fixed, i.e., state-independent tensors 
K in V: 

N^p,^j =H(C,Q,K). (3.16) 
For hyperelastic bodies with the right-invariant potential energy, 

V(<p) = V(<pB), BeO(U,r)), (3.17) 

the following holds: 

V(<p)=W(C[<p],K), (3.18) 

where again K denotes any system of state-independent tensors in V . 

An important question appears as to when the dynamics of an affine hyper- 
elastic body is simultaneously isotropic in space and matter. Obviously, this 
holds only when ()3.12j) . (|3.15| . (|3. 18|) are simultaneously satisfied. Therefore, 
the inertial tensor is spherical, J = Irj, and the potential energy V depends on tp 
only through the deformation invariants, e.g., through the quantities Ik (|1.30|) . 
or any of alternative expressions like A Q , q a , Q a or other used in the two-polar 
decomposition like ()1 .35[) . Therefore, 

V(p)=F(I 1 ,...,I n ) = G(\ 1 ,...,\ n ), (3.19) 

where, obviously, G is invariant under the group of all permutations of its 
arguments. 

It is clear that according to the general rules of Hamiltonian mechanics, in 
the potential motion of the affinely-rigid body the reactions-free affine moment 
of forces is given by 

■ dV dV u - 

dtp* A dtp A 
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and similar formulas hold for the co-moving representation. 

The relationships (|2.24p - (|2.26[) imply that in a general, not necessarily hy- 
perelastic, case equations of internal motion are simultaneously spatially and 
materially isotropic, when J — It] and N is given by (|3.9p with K = r\ or, 
equivalently, by (J3.16P with K — g. For example, in a rather academic elastic, 
but not necessarily hyperelastic, situation using the Cayley-Hamilton theorem 
one can show that 



- A 



A 

B =J2 B a(h,-.-,In)(G a -^ B , (3.21) 



where N e i b, G b are components of Net, G with the ?y-lowered index B: 

Net A b = N el AC VCB , G A B = G ag t)cb, (3.22) 

and the scalar coefficients B a in expansion (|3.21[) depend on deformation in- 
variants. One can show that in the hyperelastic case, when the potential 
V ■ ■ ■ , I n ) does exist, the coefficients are given by the following derivatives: 

dV 

B a = -2a—. (3.23) 

Obviously, the physical utility of elastic but not hyperelastic models is rather 
doubtful, nevertheless it must be admitted for the completeness of the theory. 

Another example of a doubly isotropic model is one concerning the isotropic 
internal friction in continuum droplet. The viscous stress tensor is given in a 
linear approximation by 

a% sc = 2ud?i + (C - ^) 9a b d ab g^, (3.24) 

where the constants v, ( are viscosity coefficients and cf J is the deformation 
rate tensor. In the case of affine body it is given by 

d*i = 1 (n« + ()■") , ()" = c>V/'. (3.25) 

Then, making use of the obvious formula (|3.10[) we obtain that 

det [g 



det [tp* A ] (v + fP) + (C - ^) ^ fe fc.9 lJ ) . 



K --- V °^det[ VAB ] v v nj 

(3.26) 

where Vq denotes the standard (Lagrangian) volume of the affine body. This 
agrees with the formula (|3.16l) with K = g. 

The total viscoelastic and doubly-isotropic moment of forces is given by 

N ii =N V + N % sc (3.27) 
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with the separate terms like p.2ip (|3.26p . 

Those were interesting and instructive examples of the doubly isotropic (spa- 
tially and materially) internal forces 7V IJ . It is also interesting to find a descrip- 
tion of more general isotropic forces, adapted to certain special parametrizations 
of the configuration space. In particular, some possibilities of partial separation 
of variables, or rather their subsystems, appear then. First of all, let us be- 
gin with the polar splitting (|1.38p , more precisely, with its first form where the 
orthogonal term U [<p\ stands on the left-hand side. As mentioned above, gyro- 
scopic kinetics is described by the co-moving angular velocity uj = U~ 1 dU/dt, 
whereas deformation (together with its orientation with respect to the body) is 
represented by the 77-symmetric and positive factor A [<p] in f| 1 . 38|) . As usual, 
the following tensors with 77-shifted indices will be employed: 

J A B = J AC VCB , N A B = N Ac ri CB , (3.28) 

G A B = v AC G CB , A KL = A K MV ML . (3.29) 

It is clear that 

Q = A- 1 QA + A- 1 ^ = A- 1 (q + ^A- 1 )a (3.30) 

and 

N ij = ^c^dN 00 = U l K U : > L A K cA L D N CD . (3.31) 
This suggests us to introduce the following quantity: 

N KL = A K C A L D N CD , i.e., N = (A ® A) N. (3.32) 

Just like N itself, N is also an element of U (g> U, however of a quite different 
nature. Namely, N AB are components of N with respect to the basis ipEA, 
A = 1, . . . , n, affinely co-moving with the body. Unlike this, the quantities 

j\B 

N are components of iV with respect to the orthonormal basis U [<p] Ea , 
A = 1, . . . , n, co-moving with the U [tp] -gyroscope of the polar decomposition of 

<P- 

After this substitution, our internal equations of motion, i.e., the second 
subsystem (|2.35l) . become as follows: 

aK t C a A aK rC aD auJ 0A K T C aA B r EM 
A CJ D T3 A cJ DA e 7. lA CJ D 77 — W 

dV dt dt 

+A k c J c dA d e Q e f u fm - N KM (3.33) 

with the convention (|3.28[| — (|3.29p concerning the 77-shift of tensor indices. Ob- 
viously, for the spatially isotropic models one is faced with some kind of partial 
separation of variables. Indeed, the spatial isotropy means that TV is indepen- 
dent on the variable U . It is a function of the state quantities A, dA/ dt, Q only. 
Roughly speaking, in the non-holonomic cD-representation it is a kind of cyclic 
state variable. Therefore, the procedure of solving equations of motion splits 
into three steps: 
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1) (|3.33p is a system of differential equations for the time dependence of 
quantities A, Q. 

2) Assuming that the previous step is done, we write the system of differential 
equations for U, 

^=UQ(t). (3.34) 

Let us observe that this system is time-dependent through the time evo- 
lution of u) (t) . 

3) When the steps 1), 2) are performed, we construct the final solution: 

ip(t) = U (t) A (t) . (3.35) 

Obviously, this is only the general scheme. For dynamically realistic models 
the steps 1), 2) as a rule, are not analytically solvable. Nevertheless, even 
this partial separation and a sequence of procedures is very helpful for the 
understanding the problem. In any case, its structure looks simpler and more 
adapted to operations. Let us only quote two examples corresponding to (|3.21|) . 
(|3T2%| . After substituting (f3~2g |) -(|3T2l? j) we find respectively that 

n 

Nei = ^B a (Ji.../ n )A 2a , (3.36) 



a=l 

-VadetA 



v 



dA A-l A-l dA \ 

—A 1 + A 1 — 
dt dt I 



dt 



(3.37) 



It is a nice feature of the both formulas that N e i depends only on A, and N V i SC 
depends only on A, dA/dt; there is no dependence on D. This independence 
is due to the fact that N V i SC describes the internal friction. To be honest, the 
linear dependence of N V i SC on the (^-symmetric part of Q is an approximation 
valid in the case of small internal velocities. In general, the higher powers of 
fi(w are admissible. 

Both expressions p.21[) . (|3.26p . therefore, also (|3.36p . (|3.37p . have an addi- 
tional interesting feature of being isotropic simultaneously in space and material. 
It is natural to ask for the optimal way of expressing this fact. As expected, 
the most natural way consists in using the two-polar representation (jl.36p and 
the related quantities (TPS) . ([Pi?)) , (fl~5T)|) . (fl~5"Tj) . (fT32"|) . Then, identifying the 
factors L, R in (I1.36P with linear mappings from R" to V and U respectively, 
and similarly identifying \, D with linear mappings from R" to R", we obtain 
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the following formulas: 



o = L[x + ^D- x -D$D- 1 )l \ 
<) = R[D- l xD + D- 1< ^--'4)R- 1 , (8.8!)) 



2 = fl(£-tf)iT\ (3.40) 

obviously D (|1.55p is a linear mapping from £/ to U. These formulas are simple 
and suggestive. The last of them, i.e., one for Q, is an infinitesimal expression 
of the obvious fact that U = LR^ 1 . 

As mentioned in (|2.23l) . for the potential systems the Legendre transforma- 
tion relates to K*j = K lm g mj . Similarly Y> A B is related to K ac Gcb- It is 
important that the second index is lowered with the help of Green deformation 
tensor, not with the help of the fixed material metric rjcB- One should not 
confuse Y> A b with 

K AB = (^ 1 ) A .(^ 1 ) B .if«. (3.41) 

Let us quote the explicit Legendre formulas for those quantities in the special 
materially isotropic case J AB = Ir\ AB . So, we have that 

T, = K = IL\DdD + D^- - D 2 x^j £~\ (3.42) 

S = tp^Ktp = IR [dD 2 + - DxD^j R~ l , (3.43) 

K = <fT x K {if- 1 ) 7 = mfd+^-D- 1 - D x Dj R-\ (3.44) 

Let us mention that the corresponding spin parts, i.e., doubled skew-symmetric 
parts of those quantities, equal respectively to 



S , = E-E T = IL[2DdD - D 2 x-xD 2 ) L-\ (3.45) 
V = E T — E = IR(2DxD - D 2 d -dD 2 ^j RT 1 , (3.46) 
S = K-K T = m(-DxD- 1 -D~ 1 xD + 2d) R- 1 . (3.47) 



Then for the doubly (spatially and materially) isotropic problems the quantities 
S, V are constants of motion. Unlike this, S, i.e., (|3.47p . is not a conserved 
quantity. 

In (|3 . 32(1 we have introduced the quantity N 6 U <8> V , the components of 
which represented N with respect to the moving orthonormal basis U[<p]Ea, 
A = 1, . . . , n. This representation enabled one to reduce equations of motion 
to the [/-independent form (I3.33p . Something similar may be done for the two- 
polar representation. Namely, there exist an obvious analogy between (|3.30j) 
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and (EOgD in that 

n = LQL~ 1 , n = RQR- 1 , ft = DOIT 1 . (3.48) 
It is clear that the matrix elements of the R n -tensors, 

h = X+^D~ x -DVD- 1 , Q = D- 1 xD + D- 1 ^--d, (3.49) 

are components of fi with respect to the orthonormal frame L a frozen into the 
Cauchy gyroscope and the components of fl with respect to the orthonormal 
frame co-moving with the Green deformation tensor. And the same representa- 
tion may be introduced for any other tensor quantity, in particular for the affine 
moment of forces TV. The mixed, contravariant-covariant representation of N is 
given by N%, where N — L~ 1 NL, i.e., 

N% = (L a , NL b ) = L a t N)L\, = N lk g kj . (3.50) 

Substituting the above equations to the doubly isotropic case of the internal 
subsystem (|2.35[) , we obtain the following equations of motion: 

d 2 D ,dy d& dD^ ~dD 

+D 2 X 2 - 2DxDd + Dd 2 D = ~N (d, $\ . (3.51) 

The dynamical double isotropy implies that TV depends only on the indicated 
variables D, dD/dt, £, but is independent of the angular variables L. R. 
Therefore, similarly like in (|3.33p . there is a partial separability of the problem 
HMD: 

1) Just as it was the case with p.33|) . one solves the system p. 511) . To 
be more precise, one dreams about solving this system of n 2 ordinary 
differential equations for the n 2 dynamical variables D, £, Some kind 
of rigorous solutions is possible only for the two-dimensional case n = 2. 
For higher dimensions, including the physical case n = 3, only some special 
solutions may be analytically found. 



2) When the time dependence M3f4 [x(t)> J ls "known", we substitute 
it to the definition of angular velocities: 

dL dR -j. . 

Tt =L ^ Tt=»- (3 - 52) 

Then one obtains the system of differential equations with right-hand sides 
explicitly dependent on time. 

3) After "solving" (|3.52p we substitute everything to (|1.36p and obtain the 
final solution. 

As mentioned many times above, this partial reduction (separability) of (|3.33|) . 
(|3.5ip is rather ideal and qualitative, nevertheless, it is helpful in understanding 
the dynamical structure of spatially and doubly isotropic models. 
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4 D'Alembert and vakonomic models of rotation- 
less motion 



Let us now discuss briefly the interesting special case of non-holonomic rotation- 
less constraints. This will be rather an introductory analysis; up to our knowl- 
edge nobody discussed this kind of constraints, either in the d'Alembert or 
vakonomic version. From some point of view the apparently exotic vakonomic 
form is rather simpler and more elegant [6]. It is yet rather too early to try 
deciding which is more physical and in what kind of problems. 

Let us substitute formally the polar representation (|3.4j) . (|3.5j) of (|2.61j) to 
the polar expression of the kinetic energy (|2.4[) . Then we obtain that 

1 dA K A dA L B AB 1 , nK dA D c AC dA L B AB 

Tint = o^KL—t: "77 + -jVKL (A ) D A A J 

8 dt dt 4 v ' dt dt 



Ivkl (A~ 1 ) K e—^—A c a (A-') L F ^A D B J AB . (4.1) 



After calculations this may be expressed in the following more concise form: 
T mt = l -mL(A^) K E {A-') L F (A E c dA ' A 



dt 

dAE "A° A ) + d 4±A% ) J AB . ,4.2, 



dt J \ dt dt 

The variational derivative of Ti nt with respect to the symmetric tensor 

A A b = tiacA c b = Aba (4.3) 



is given by 



SA AB 



sy mm 



1 d 

'dt V dt 



riKL^\A^(A-^ A A B ) D ]j ED 



)mA ( (A-'f^A^ (A-y^ D ) J FD 



4 1 dt y ' dt 
1 dA K E dA F 
A VKL — t dt 



1 dA K E dA F D D (a -is L (AaB) t eg 
iIkl—t, j— A G {A ) A > F J 



\vkl (A-Y E ^A C M ^A% (A-Y (A A B W MN 

K dAK °(A-^ L dAE(A JB)D 

1 / ._ n jf dA E c c 



1 , , ,,k dA E c .n , . ni dA F( - A mn , . 
-;Vkl(A- 1 ) e —^A c d (A-i) f——J b)d . (4.4) 
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When there are hyperelastic forces derivable from the potential V depend- 
ing only on the Green deformation tensor, then equations of motion have the 
following form: 

ST "" = -A KC v KiA N B ^ c , (4.5) 

symm 



SA AB 

where 

n bu = _ ( DV y^ (4, 6 ) 

In spite of their apparently complicated structure, equations (|4.4| are read- 
able. And having them solved for the time dependence of Aab, we obtain from 
(13. 4p / (13. 5p the time dependence of a), and then, solving (in principle) (|1.55|) for 
dependence t — > U(t), we finally obtain (in principle) ip = UA. 

Let us mention that all tensor indices are shifted from their natural position 
with the help of rj. 

The usual d'Alembert procedure, i.e., the symmetric part of (|3.33|) with al- 
gebraically substituted constraints (|2.61[) . i.e., (|3.4p / (|3.5p , leads to the following 
form, less readable than ([P|). (14~51) . (|4"1?)) : 



J ~d^ A )a ~ JbA E J t 2 V A > F Jt\ A ) A 



D) A 



(A e f) (A'Y iC A D \) - J A b^ ( {A-') E A (a«°) A d \ 



dt J dt \ dt 

-iiA'eHA- 1 ) *) {^) F 4A AH(C ) AD)A 

-(AKii^A) = N (CD) . (4.7) 

Again the 77-shift of indices is meant here. The difference between (|4.4p / (14. 5[) 
and (|4.7p on their right-hand side is not essential, because it is only due to the 
A-term transformation of N into N . They may be written in a similar form in 
this sense. But the difference between other terms of (|4.6|) and (|4.7I1 is more 
essential. The detailed analysis of this difference is postponed to the next paper. 
In any case, it is a general rule that the d'Alembert and vaconomic procedures 
give different equations. 



Appendix: d'Alembert vs. vaconomic constraints 

The problem appear more than century ago. It is well known that when the 
holonomic constraints 

F a (q)=0, o = l,...,m, (4.8) 

are imposed onto the motion of a Lagrangian dynamical system with general- 
ized coordinates q , . . .,5™, then one can equivalently use the d'Alembert pro- 
cedure or the restricted extremum (more precisely, stationary value) problem. 
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If Lagrangian is given by L(q,q), then the Lusternik theorem tells us that the 
conditional extremum (more precisely, stationary value) 

sJldt^O, F a (q) = 0, a=l,...,m, (4.9) 

is given by the functions of time satisfying equations: 



where 

And those are exactly d'Alembert equations with the multipliers 

A" . ,4.12) 

The same formulas for reaction forces Ri hold also for non- variational, e.g., 
dissipative dynamical models: 

d dL dL n „ o , 

jtw-w =Q - (413) 

where Qi are non- variational generalized forces. Then as well we have that 

So, there was a natural temptation to expect something similar for systems with 
non-holonomic constraints, for simplicity linear in velocities, 

w -i(9)^=«, (4-15) 
but without the intergrability assumption 

OF 

"ai^ = 0, (4.16) 
i.e., without the vanishing of exterior differentials: 

But it turned out in contrary: d'Alembert procedure gives again the equations 

d dL dL . .dq % , , H „, 

- w - w = Q i + X- Uait uM ^=0 t (4.18) 
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with reactions coefficients A a to be eliminated. But the Lusternik theorem for 

&jLdt = {), u ai (q)^=0 (4.19) 

gives something drastically else: 

d dL dL d\i a f duj ai dui a j\ .,■ n on s 

dd 1 

= 0. (4.21) 

The difference is obvious. Moreover, fi a become a kind of dynamical variables, 
because they occur both by itself and their time derivatives. The problems 
of sliding-free rolling motion are ruled by the d'Alembert procedure. But, on 
the other hand, the Lusternik variational, i.e., vaconomic, procedure looks very 
interesting and intriguing. It gives rice to the new mathematical discipline and 
its applications seem to be also possible, first of all, in active control problems. 

Our equations for the rolling-free affine motion in the d'Alembert and va- 
conomic sense are also drastically different, although as yet we are unable to 
express them in qualitative terms. 
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